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Biodiversity is essential to the viability of ecological systems. Species diver- 
sity in ecosystems is promoted by cyclic, non-hierarchical interactions among 
competing populations. Central features of such non-transitive relations are 
represented by the 'rock-paper-scissors' game, where rock crushes scissors, scis- 
sors cut paper, and paper wraps rock. In combination with spatial dispersal of 
static populations, this type of competition results in the stable coexistence of 
all species and the long-term maintenance of biodiversity^^. However, popula- 
tion mobility is a central feature of real ecosystems: animals migrate, bacteria 
run and tumble. Here, we observe a critical influence of mobility on species 
diversity. When mobility exceeds a certain value, biodiversity is jeopardized 
and lost. In contrast, below this critical threshold all subpopulations coexist 
and an entanglement of travelling spiral waves forms in the course of time. We 
establish that this phenomenon is robust, it does not depend on the details of 
cyclic competition or spatial environment. These findings have important impli- 
cations for maintenance and temporal development of ecological systems and are 
relevant for the formation and propagation of patterns in microbial populations 
or excitable media. 

The remarkable biodiversity present in ecosystems confounds a naive interpretation of 
Darwinian evolution in which interacting species compete for limited resources until only 
the fitter species survives. As a striking example, consider that a 30-g sample of soil from 
a Norwegian forest is estimated to contain some 20,000 common bacterial species^. Evo- 
lutionary game theory^"^, in which the success of one species relies on the behaviour of 
others, provides a useful framework to investigate co-development of populations theoreti- 
cally. In this context, the rock-paper-scissors game has emerged as a paradigm to describe 
species If three subpopulations interact in this non-hierarchical way, one 

expects intuitively that diversity may be preserved: Each species dominates another only 
to be outperformed by the remaining one in an endlessly spinning wheel of species-chasing- 
species. 

Communities of subpopulations exhibiting such dynamics have been identified in numer- 
ous ecosystems, ranging from coral reef invertebrates^^ to lizards in the inner Coast Range of 



California^. In particular, recent experimental studies using microbial laboratory cultures 
have been devoted to the influence of spatial structure on time development and coexistence 
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of species^^-H Investigating three strains of colicinogenic E. coli in different environments, it 



has been shown that cyclic dominance alone is not sufficient to preserve biodiversity. Only 
when the interactions between individuals are local (e.g. bacteria arranged on a Petri dish), 
spatially separated domains dominated by one subpopulation can form and lead to stable 



coexistence^]. 

In this Letter, we show that biodiversity is affected drastically by spatial migration of indi- 
viduals, a ubiquitous feature of real ecosystems. Migration competes with local interactions 
such as reproduction and selection, thereby mediating species preservation and biodiversity. 
For low values of mobility, the temporal development is dominated by interactions among 
neighbouring individuals, resulting in the long-term maintenance of species diversity. In 
contrast, when species mobility is high, spatial homogeneity results and biodiversity is lost. 
Interestingly, a critical value of mobility sharply delineates these two scenarios. We obtain 
concise predictions for the fate of the ecological system as a function of species mobility, 
thereby gaining a comprehensive understanding of its biodiversity. 

The influence of mobility on species coexistence was previously studied within the frame- 
work of coupled habitat patches ("island models" jlSHISl j n particular, Levin considered 
an idealized two-patch system and observed a critical mobility for stable coexistence^^. 
Other models comprising many spatially arranged patches were shown to facilitate pattern 
formatioE pP^ . As often in nature spatial degrees of freedom vary continuously (e.g. bacte- 
ria can visit the entire area of a Petri dish), we relax the simplifying assumption of habitat 
patches and consider continuous spatial distribution of individuals. Moreover, as an inherent 
feature of real ecosystems and in contrast to previous deterministic invest igat ion JEMUS we 
explicitly take the stochastic character of the interactions among the populations into ac- 
count. Such interacting particle systems, where individuals are discrete and space is treated 
explicitly, were already considered in ecological contexts - 1 2 | 4 | 5 | 20 ] rj ne ^g^yjoyj. Q f these 



models often differs from what is inferred from deterministic reaction-diffusion equations, 



or from interconnected patches^. In the case of cyclic competition, such stochastic spatial 
systems have been shown to allow for stable coexistence of all species HZEE] w hen individuals 
are static. Here we explore the novel features emerging from individuals' mobility. 

Consider mobile individuals of three subpopulations (referred to as A, B, and C), ar- 
ranged on a spatial lattice, where they can only interact with nearest neighbours. For 
the possible interactions, we consider a version of the rock-paper-scissors game, namely 
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a stochastic spatial variant of the model introduced in 1975 by May and Leonard^ ( see 
Methods). Schematic illustrations of the model's dynamics are provided in Fig. 1. The 
basic reactions comprise selection and reproduction processes, which occur at rates a and 
/i, respectively. Individuals' mobility stems from the possibility for two neighbouring indi- 
viduals to swap their position (at rate e) and to move to an adjacent empty site. Thereby, 
individuals randomly migrate on the lattice. We define the length of the square lattice as 
the size unit, and denote by N the number of sites. Within this setting, and applying the 
theory of random walkJ^, the typical area explored by one mobile individual per unit time 
is proportional to M = 2eiV _1 , which we refer to as the mobility. The interplay of the latter 
with selection and reproduction processes sensitively determines whether species can coexist 
on the lattice or not, as discussed below. 

We performed extensive computer simulations of the stochastic system (see Methods) 
and typical snapshots of the steady states are reported in Fig. 2. When the mobility of the 
individuals is low, we find that all species coexist and self-arrange by forming patterns of 
moving spirals. Increasing the mobility M, these structures grow in size, and disappear for 
large enough M. In the absence of spirals, the system adopts a uniform state where only 
one species is present, while the others have died out. Which species remains is subject to 
a random process, all species having equal chances to survive in our model. 

Concise predictions on the stability of three-species coexistence are obtained by adapting 
the concept of extensivity from statistical physics (see Supplementary Information). Namely, 
we consider the typical waiting time T until extinction occurs, and its dependence on the 
system size N. If T(N) ~ N, the stability of coexistence is marginal^. Conversely, longer 
(shorter) waiting times scaling with higher (lower) powers of N indicate stable (unstable) 
coexistence. These three scenarios can be distinguished by computing the probability P cxt 
that two species have gone extinct after a waiting time t ~ N. In Fig. 2, we report the 
dependence of P ext on the mobility M. For illustration, we have considered equal reaction 
rates for selection and reproduction, and, without loss of generality, set the time-unit by 
fixing a — \i — 1. Increasing the system size N, a sharpened transition emerges at a critical 
value M c = (4.5 ± 0.5) x 10~ 4 for the fraction of the entire lattice area explored by an 
individual in one time-unit. Below M c , the extinction probability P cxt tends to zero as the 
system size increases, and coexistence is stable (implying super-persistent transients^, see 
Supplementary Information). On the other hand, above the critical mobility, the extinction 
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probability approaches one for large system size, and coexistence is unstable. As a central 
result of this Letter, we have identified a mobility threshold for biodiversity: 

There exists a critical value M c such that a low mobility M < M c guarantees coexistence 
of all three species, while M > M c induces extinction of two of them, leaving a uniform state 
with only one species. 

To give a biological illustration of this statement, let us consider colicinogenic strains of 
E.coli growing on a Petri disfM In this setting, 10 bacterial generations have been observed 
in 24 hours, yielding selection and reproduction rates of about 10 per day. As the typical 
size of a Petri dish is roughly 10 cm, we have evaluated the critical mobility to be about 
5 x 10 2 /im 2 /s. Comparing that estimate to the mobility of E.coli, we find that it can, by 
swimming and tumbling in super soft agar, explore areas of more than 10 3 /im 2 per second^. 
This value can be considerably lowered by increasing the agar concentration. 

When the mobility is low, i.e. M < M c , the interacting subpopulations exhibit fascinating 
patterns, as illustrated by the snapshots of Fig. 2. The emerging reactive states are formed 
by an entanglement of spiral waves, characterising the competition among the species which 
endlessly hunt each other, as illustrated in movies 1 and 2 (see Supplementary Informa- 
tion). Formation of this type of patterns has been observed in microbial populations, such 
as myxobacteria aggregation^ or multicellular Dictyostelium mounds^, as well as in cell 
signaling and control^. Remarkably, a mathematical description and techniques borrowed 
from the theory of stochastic processed^ allow to obtain these complex structures by means 
of stochastic partial differential equations (PDE), see Fig. 3 and Methods. Furthermore, 
recasting the dynamics in the form of a complex Ginzburg-Landau equatior P^"^ allows to 
obtain analytical expressions for the spirals' wavelength A and frequency (see Supplemen- 
tary Information). These results, up to a constant prefactor, agree with those of numerical 
computations; see our forthcoming publication (T.R., M.M. and E.F., in preparation). 

As shown in Fig. 2, the spirals' wavelength A rises with the individuals' mobility. Our 
analysis reveals that the wavelength is proportional to \f~M (see Supplementary Informa- 
tion). This relation holds up to the mobility M c , where a critical wavelength A c is reached. 
For mobilities above the threshold M c , the spirals' wavelength A exceeds the critical value 
A c and the patterns outgrow the system size causing the loss of biodiversity, see Fig. 2. We 
have found A c to be universal, i.e. independent on the selection and reproduction rates. 
This is not the case for M c , whose value varies with these parameters (see Supplementary 
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Information). Using lattice simulations, stochastic PDE and the properties of the complex 
Ginzburg-Landau equation, we have derived the dependence of the critical mobility M c (/x) 
on the reproduction rate /i (where the time-unit is set by keeping a = 1). It enables to 
analytically predict, for all values of parameters, whether biodiversity is maintained or lost. 
We have summarized these results in a phase diagram, reported in Fig. 4. One identifies a 
uniform phase, in which two species go extinct (when M > M c (/i)), and a biodiverse phase 
(when M < M c (/x)) with coexistence of all species and propagation of spiral waves. 

The generic ingredients for the above scenario to hold are the mobility of the individuals 
and a cyclic dynamics exhibiting an unstable reactive fixed point. The underlying mathemat- 
ical description of this class of dynamical systems is in terms of complex Ginzburg-Landau 
equations. Their universality classes reveal the robustness of the phenomena which we have 
reported above, i.e. the existence of a critical mobility and the emergence of spiral waves; 
they are not restricted to specific details of the model. 

Our study has direct implications for experimental research on biodiversity and pattern 
formation. As an example, one can envisage an experiment extending the study^ on coli- 
cinogenic E.coli. Allowing the bacteria to migrate in soft agar on a Petri dish should, for 
low mobilities, result in stable coexistence promoted by the formation of spiral patterns. 
Increasing the mobility (e.g. on super soft agar), the patterns should grow in size and finally 
outgrow the system at some critical value, corresponding to the threshold M c discussed 
in this Letter. For even higher values of the mobility, biodiversity should be lost after a 
short transient time and only one species should cover the entire Petri dish. We think that 
both mobility's regimes, corresponding to the biodiverse and uniform phases, should be 
experimentally accessible. 

We have shown how concepts from game theory combined with methods used to study 
pattern formation reveal the subtle influence of mobility on the temporal development of 
coexisting species. Many more questions and applications regarding the seminal interplay 
between these different fields lie ahead. As an example, concerning the evolution of coop- 
eration, it has been shown that cyclic dominance can occur in social dilemmas 9 I 3CH which 
suggests implications of our results to behavioural sciences. 
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Methods 



To model cyclic dominance, we use a stochastic lattice version (following work by Durrett 
and Levin) of a model proposed by May and Leonard in 1975. As main characteristics, in 
the absence of spatial structure, their equations possess a deterministically unstable fixed 
point associated to coexistence of all three species: In the course of time, the system spirals 
(in the phase space) away from coexistence and moves in turn from a state with nearly only 
A's to another one with nearly only B's, and then to a state with almost only C's. 

In our stochastic lattice simulations, we have arranged the three subpopulations on a two- 
dimensional square lattice with periodic boundary conditions. Every lattice site is occupied 
by an individual of species A, B, or C, or left empty. At each simulation step, a random 
individual is chosen to interact with one of its four nearest neighbours, which is also randomly 
determined. Whether selection, reproduction or mobility occurs, as well as the corresponding 
waiting time, is computed according to the reaction rates using an efficient algorithm due to 
Gillespie^l We set one generation, i.e. when every individual has reacted on average once, 
as unit of time. To compute the extinction probability, we have used different system sizes, 
from 20 x 20 to 200 x 200 lattice sites, and sampled between 500 and 2000 realizations. The 
snapshots shown in Fig. 2 result from system sizes of up to 1000 x 1000 sites. 

Our stochastic PDE consist of a mobility term, nonlinear terms describing the determinis- 
tic temporal development of the nonspatial model (May-Leonard equations), and (multiplica- 
tive) white noise, see Supplementary Information. We have solved the resulting equations 
with the help of open software from the XMDS projectP^, using the semi-implicit method 
in the interaction picture (SUP) as an algorithm, spatial meshes of 200 x 200 to 500 x 500 
points, and 10, 000 points in the time direction. 

Supplementary Information 

is available online. 
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Figure captions 



Figure 1 

The rules of the stochastic model. Individuals of three competing species A (red), B (blue), and 
C (yellow) occupy the sites of a lattice, a, They interact with their nearest neighbours through 
selection (i) or reproduction (ii). Both reactions occur as Poisson processes at rate a and fj,, 
respectively. Selection reflects cyclic dominance: A can kill B, yielding an empty site (black) 
there. In the same way, B invades C, and C in turn outcompetes A. Reproduction of individuals 
is only allowed on empty neighbouring sites, mimicking a finite carrying capacity of the system. We 
also endow individuals with mobility: at rate e, they are able to swap position with a neighbouring 
individual or hop on an empty neighbouring site (iii). b, An example of the processes (i)-(iii), 
taking place on a 3 x 3 square lattice. 

Figure 2 

The critical mobility M c . Mobility below the value M c induces biodiversity; while it is lost above 
that threshold, a, We show snapshots obtained from lattice simulations of typical states of the 
system after long temporal development (i.e. at time t ~ N) and for different values of M (each 
color, blue, yellow and red, represents one of the species and black dots indicate empty spots). 
Increasing M (from left to right), the spiral structures grow, and outgrow the system size at the 
critical mobility M c : then, coexistence of all three species is lost and uniform populations remain 
(right), b, Quantitatively, we have considered the extinction probability P ex t that, starting with 
randomly distributed individuals on a square lattice, the system has reached an absorbing state 
after a waiting time t = N. We compute P cx t as function of the mobility M (and a = fi = 1), 
and show results for different system sizes: N = 20 x 20 (green), N = 30 x 30 (red), N = 40 x 40 
(purple), N = 100 x 100 (blue), and N = 200 x 200 (black) . As the system size increases, the 
transition from stable coexistence (Pext = 0) to extinction (P cx t = 1) sharpens at a critical mobility 
M c w (4.5 ±0.5) x 10~ 4 . 
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Figure 3 

Spiralling patterns, a, Typical spiral (schematic). It rotates around the origin (white dot) at a 
frequency uo and possesses a wavelength A. b, In our lattice simulations, when the mobility of 
individuals lies below the critical value, all three species coexist, forming mosaics of entangled, 
rotating spirals (each color represents one of the species and black dots indicate empty spots), 
c, We have found that the system's development can aptly be described by stochastic partial 
differential equations. In the case of lattice simulations and stochastic partial differential equations, 
internal noise acts as a source of local inhomogeneities and ensures the robustness of the dynamical 
behaviour: the spatio-temporal patterns are independent of the initial conditions, d, Ignoring the 
effects of noise, one is left with deterministic partial differential equations which also give rise 
to spiralling structures. The latter share the same wavelength and frequency with those of the 
stochastic description but, in the absence of fluctuations, their overall size and number depend on 
the initial conditions and can deviate significantly from their stochastic counterparts. In b and c, 
the system is initially in a homogeneous state, while d has been generated by considering an initial 
local perturbation. Parameters are a = \i = 1 and M = 1 X 10~ 5 . 

Figure 4 

Phase diagram. The critical mobility M c as a function of the reproduction rate [i yields a phase 
diagram with a phase where biodiversity is maintained as well as a uniform one where two species 
go extinct. Time unit is set by a = 1. On the one hand, we have computed M c from lattice 
simulations, using different system sizes. The results are shown as blue crosses. On the other 
hand, we have calculated M c using the approach of stochastic PDE (black dots, black lines are 
a guide to the eye) as well as analytically via the complex Ginzburg-Landau equation (red line). 
Varying the reproduction rate, two different regimes emerge. If (i is much smaller than the selection 
rate, i.e. [i <C <r, reproduction is the dominant limiter of the temporal development. In this case, 
there is a linear relation with the critical mobility, i.e. M c ~ /i, as follows from dimensional analysis. 
In the opposite case, if reproduction occurs much faster than selection (|j > it), the latter limits 
the dynamics and M c depends linearly on a, i.e. M c ~ a. Here, as a = 1 is kept fixed (time-scale 
unit), this behaviour reflects in the fact that M c approaches a constant value for fj, S> a. 
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a (i) Selection, rate a: (ii) Reproduction, rate \i: 
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(ii) Reproduction (fj.) 




(iii) Exchange (e) 



14 



Figure 2 




15 



Figure 3 



a Typical spiral b Lattice simulations c Stochastic PDE d Deterministic PDE 




1(3 




17 



Mobility promotes and jeopardizes biodiversity in 
rock-paper-scissors games 

Tobias Reichenbach, Mauro Mobilia, and Erwin Frey 
Supplementary Information 

In this Supplementary Information, we further elaborate our analysis by explaining some 
technical aspects of the Letter and illustrate our findings by providing two supplementary 
movies. The latter illustrate the spatio-temporal dynamics of the biodiverse reactive states 
occurring in the system under consideration. 

Below, we first present the concept of extensivity, which allows to precisely discriminate 
between the regime where biodiversity is stable (i.e. maintained) and the situation where it 
is unstable and the system settles in one of the uniform (absorbing) states. Then, we present 
details on the stochastic differential equations as well as the complex Ginzburg-Landau 
equation used to analyse the system. 

We also show how the spirals' wavelength A is related to the mobility M. We explain 
that such a relation has allowed to derive the functional dependence M c (/x) of the critical 
mobility on the reproduction rate. Finally, the main findings reported in the Letter are 
revisited in a supplementary discussion centered on the information conveyed by the movies. 

Extensivity 

Even if coexistence appears stable, as observed for low mobilities, there is a certain 
probability that two species go extinct due to possible large yet rare fluctuations. Indeed, 
the only absorbing states where no reactions (and therefore no fluctuations) occur, are the 
uniform configurations where only one species survives. For this reason, these are the only 
stable states in the long run. However, the typical waiting time T(N) until extinction 
occurs is generally very long when the system size N is large. This suggests to consider 
the dependence of the waiting time T(N) on N. Quantitatively, we discriminate between 
stable and unstable coexistence by using the concept of extensivity, adapted from statistical 
physics. If the ratio T/N tends to infinity (T/N — > oo) in the asymptotic limit of large 
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systems (N — > oo), the typical waiting time strongly prolongs with N (typically with an 
exponential dependence). This scenario is called super- extensive or stable. On the other 
hand, if T/N — > 0(1) (i.e. the ratio approaches a finite non-zero value) that is referred 
to as the extensive case, which has been shown to correspond to marginal (or neutral) 
stability^. Instability of the coexistence state (towards a uniform one) is encountered 
when T/N — ► (sub-extensive scenario), where the waiting time is short as compared to 
the system size. 

These definitions of stability and instability (with neutral stability separating the two cases) 
in the presence of absorbing states are intimately related to the concept of transients. In 
fact, HastingJ^ suggested to study not only the ultimate fate of a system, but also to 
consider the behaviour at smaller (and probably ecologically more relevant) time scales. 
According to the definition introduced above, stable or neutrally stable coexistence implies 
coevolution of the populations for a very large number of generations. This corresponds to 
the existence of extremely long-lived persistent transients^ (super-persistent). It is also 
worth noticing that transients lasting for several generations can occur even in the case of 
unstable coexistence. This typically happens when the number of individuals N is large. 
In the situation of Fig. 2, we have considered the extinction probability P cxt that, starting 
from random initial conditions (i.e. spatially homogeneous configurations, with equal 
concentrations of each species), the system has reached a uniform state after a time t 
proportional to the system size, i.e. t ~ N. In the asymptotic limit N — > oo, three distinct 
cases arise. In a first regime, the extinction probability tends to zero with the system size 
N. In Fig. 2, this occurs when M < M c . This scenario corresponds to the superextensive 
situation (i.e. T/N — > oo, with N — > oo) where the coexistence of all populations is stable. 
As a second case, the extinction probability approaches a finite value between and 1, i.e. 
T/N — > 0(1), and we recover neutral stability. In Fig. 2, such a behaviour arises exclusively 
at the vicinity of the critical mobility M c . In a third regime, the extinction probability does 
reach the value 1, which means that T/N — > 0. This is the subextensive scenario where the 
coexistence is unstable and biodiversity is lost. In Fig. 2, this happens above the critical 
mobility, i.e. when M > M c . 
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Stochastic partial differential equations 



Within the theory of stochastic processes^, the dynamics of the stochastic lattice system 
is described by a master equation. In the limit of large systems, using a Kramers-Moyal 
expansion, the latter allows for the derivation of a proper Fokker-Planck equation, which in 
turn is equivalent to a set of stochastic partial differential equations. The latter consist of 
a mobility term, nonlinear terms describing the deterministic dynamics of the nonspatial 
model (May-Leonard equations), and noise terms. For the noise terms, we have found 
that contributions stemming from selection and reproduction events scale as iV -1 / 2 , while 
fluctuations originating from exchanges (mobility) decay as A r_1 ; the latter may therefore 
be ignored. What remains, is (multiplicative) white noise whose strength scales as N~ l l 2 . 

Complex Ginzburg-Landau Equation (CGLE) 

Ignoring the noise terms in the stochastic differential equations describing the system, 
the resulting partial differential equations fall into the class of the Poincare-Andronov-Hopf 
bifurcation, known from the mathematics literature^. Applying the theory of center mani- 
folds and normal forms developed there, we have been able to cast the deterministic equations 
into the form of the complex Ginzburg-Landau equation: 



where z is a complex variable and c±, C3 are constants depending on the rates a and /1. This 
equation leads to the formation of dynamic spirals and allows to derive analytic results for 



Scaling relation and critical mobility 

An important question is to understand what is the mechanism driving the transition 
from a stable coexistence to extinction at the critical mobility M c . To address this issue, 
we first note that varying the mobility induces a scaling effect, as illustrated in Fig. 2. In 
fact, increasing the mobility rate M results in zooming into the system. As discussed above 
(see the main text and Methods), the system's dynamics is described by a set of suitable 




(1) 



their wavelength and frequency, see e.g. the review by Aranson and Kramer^! 
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stochastic partial differential equations (SPDE) (T.R., M.M., and E.F., in preparation) 
whose basic properties help rationalize this scaling relation. In fact, the mobility enters the 
stochastic equations through a diffusive term MA, where A is the Laplace operator involving 
second-order spatial derivatives. Such a term is left invariant when M is multiplied by a 
factor a while the spatial coordinates are rescaled by y/a. It follows from this reasoning 
that varying M into ctM translates in a magnification of the system's characteristic size by 
a factor y/a (say a > 1). This implies that the spirals' wavelength A is proportional to \/M 
(i.e. A ~ y/~M) up to the critical M c . 

When the spirals have a critical wavelength A c , associated with the mobility M c , these 
rotating patterns outgrow the system size which results in the loss of biodiversity (see the 
main text). In the "natural units" (length is measured in lattice size units and the time-scale 
is set by keeping a — 1), we have numerically computed A c = 0.8 ± 0.05. This quantity has 
been found to be universal, i.e. its value remains constant upon varying the rates a and fi. 
However, this is not the case of the critical mobility M c , which depends on the parameters of 
the system. Below the critical threshold M c , the dynamics is characterized by the formation 
of spirals of wavelength A(/i, M) ~ \/M. This relation, together with the universal character 
of A c , leads to the following equation: 

m ^ = (auOT)) 2m ' < 2) 

which gives the functional dependence of the critical mobility upon the system's parameter. 
To obtain the phase diagram reported in Fig. 4 we have used Eq. (J2J) together with values 
of A(/i, M) obtained from numerical simulations. For computational convenience, we have 
measured A(/i, M) by carrying out a careful analysis of the SPDE's solutions. The results 
are reported as black dots in Fig. 4. We have also confirmed these results through lattice 
simulations for systems with different sizes and the results are shown as blue dots in Fig. 2. 
Finally, we have taken advantage of the analytical expression (up to a constant prefactor, 
taken into account in Fig. 2) of A(/x, M) derived from the complex Ginzburg-Landau equa- 
tion (CGLE) associated with the system's dynamics (see Methods): with Eq. (|2j), we have 
obtained the red curve displayed in Fig. 2. This figure corroborates the validity of the var- 
ious approaches (SPDE, lattice simulations and CGLE), which all lead to the same phase 
diagram where the biodiverse and the uniform phases are identified. 
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Supplementary Movie 1 

In the first movie, the dynamics of individuals of species A, B and C follows the reactions 
illustrated in Fig. 1 with rates \i — 1 (reproduction), a = 1 (selection) and e = 2.4 (exchange 
rate). In Movie 1, individuals of each species are indicated in different colours (empty sites 
are shown as black dots). The dynamics takes place on a square lattice of N = 400 x 400 
sites, such that there are up to 1.6 x 10 5 individuals in the system. This set of parameters 
corresponds to a mobility rate M = 2e/N = 3 x 10~ 5 well below the critical threshold 
M c f=s 4.5 ± 0.5 x 10~ 4 (see Figs. 2, 4 and text). Initially the system is in a well-mixed 
configuration with equal density of individuals of each species and empty sites. As time 
increases and since M < M c , biodiversity is maintained and complex dynamical patterns 
form in the course of the temporal development resulting in a rich entanglement of spiral 
waves. 

Supplementary Movie 2 

In the second movie, the mobility of the individuals has been increased. In fact, the 
dynamics of individuals of species still follows the reactions illustrated in Fig. 1 with rates 
/i = 1 (reproduction) and a = 1 (selection), but the exchange rate is now e = 6. In Movie 2, 
individuals of each species are still indicated in different colours (empty sites are shown as 
black dots). The dynamics takes place on a square lattice of N = 200 x 200 sites, allowing up 
to 4 x 10 4 individuals in the system. This set of parameters corresponds to a mobility rate 
M = 3 x 10~ 4 relatively close to the critical threshold M c « 4.5 ± 0.5 x 10~ 4 (see Figs. 2, 4 
and text). Initially the system is in a well- mixed state with equal density of individuals 
of each species and empty sites. As time increases and since M < M c , biodiversity is still 
maintained and patterns form in the course of the time development. However, as compared 
to the first movie, one notices that the size of the patterns has increased and one now only 
distinguishes one pair of antirotating spirals. 

Supplementary Discussion 
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The supplementary movies illustrate the system's time development in the coexistence 
phase, i.e. the emergence of dynamical complex patterns deep in that phase (Movie 1) and 
close to (yet below) the threshold M c (Movie 2, see text and Fig. 3). 

Starting from initially homogeneous (well-mixed) configurations, after a short transient 
regime, spiral waves rapidly emerge and characterize the long-time behaviour of the sys- 
tem which settles in a reactive steady state (super- extensive case, see text). The wavefronts, 
merging to form entanglement of spirals, propagate with spreading speed v* and wavelength 
A. In Movies 1 and 2, it appears clearly that by rising the individuals' mobility, one increases 
the wavefronts propagation velocity and the wavelength of the resulting dynamical patterns, 
as well as the size of each spiral. From PDE associated with the system's dynamics, we can 
rationalize this discussion and estimate these quantities for the cases illustrated in Movies 
1 and 2. Namely, for the spreading speed, we have obtained v* ~ 3.5 x 10 -3 (lattice-size 
units per time-step, Movie 1) and v* ~ 1.1 X 10~ 2 (Movie 2). Similarly, the wavelength were 
found to be A ~ 0.21 (lattice-size units, Movie 1) and A ~ 0.65 (Movie 2). Here, rising the 
rate M from 3 x 10 -5 (Movie 1) to 3 x 10~ 4 (Movie 2) results in the enhancement of v* 
and A by a factor V^IO ~ 3.16. In Movie 2, the size of the spirals can also be estimated to 
have been magnified by the same factor -\/T0 ~ 3.16 with respect to those of Movie 1. As 
explained in the text, this scaling property of the system can be understood by considering 
the stochastic partial differential equation describing the dynamics, which were obtained 
from the underlying master equation through a system size expansion (see Methods). 
By rising the individuals' mobility, one therefore increases the size of the spiralling patterns 
(whose wavelength is proportional to \fM) and for sufficiently large value of the exchange 
rate (i.e. of M), as in Movie 2, just a few spirals nearly cover the entire lattice. This happens 
up to the critical value A c ~ 0.8, found to be universal. In fact, when A > A c the whole 
system is covered with one single ("giant") spiral which cannot fit within the lattice. This 
effectively results in the extinction of two species and the loss of biodiversity. As explained 
in the text, by exploiting the fact that A cx vM and the universal character of A c , one can 
infer the existence of the critical mobility rate M c = M c (/i) [see Eq. (J2}], as illustrated in 
Fig. 4. This allows to discuss the fate of the system (i.e. biodiversity versus extinction) in 
terms of the reaction and mobility rates fi and M, respectively: For given reaction rates /i 
and e (without loss of generality a is set to unity, see text) and system size L, we obtain 
a critical value M c (/i) of the mobility rate. In fact, a reactive steady state is reached (and 
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biodiversity maintained) only if M < M c (fi). When the individuals' mobility is too fast, 
i.e. when M > M c (/x), the system can be considered to be well-mixed and its dynamics 
therefore can be aptly described in terms of homogeneous rate equations which predicts the 
extinction of two species (see Methods). 

In the cases illustrated in Movies 1 and 2, M c m 4.5 ± 0.5 x 10~ 4 and the wavefronts prop- 
agate with A < A c , so that biodiversity is always preserved. However, we notice that the 
resulting spatio-temporal patterns are quite different: while one finds a rich entanglement 
of spirals deep in the coexistence phase (i.e. for M pa 3 x 10~ 5 <C M c , Movie 1), only one 
pair of antirotating spirals fill the system when one approaches the critical value M c (Movie 
2). 
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